#!/bin/bash -e

function info() {
echo Usage: `basename $0` hom_bulk.txt het_single.txt hom_single.txt
exit 65
}

while getopts  ":i:b:p:" opts
do
        case  $opts  in
        i) interval=$OPTARG;;
		b) asemmbly=$OPTARG;;
		p) out_prefix=$OPTARG;;
		*) info;;
        esac
done
shift $(($OPTIND - 1))

if [ -z $1 ]; then info; fi


. /mnt/ilustre/app/medical/tools/.var

Rscript $tools_path/script/fpr.r $@

rename 1.fpr.txt $out_prefix.fpr.txt 1.fpr.txt